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ABSTRACT 

A deep search for the potential glycine precursor hydroxylamine (NH 2 OH) using the 
Caltech Submillimeter Observatory (CSO) at A = 1.3 mm and the Combined Array for 
Research in Millimeter-wave Astronomy (CARMA) at A = 3 mm is presented toward 
the molecular outflow L1157, targeting the B1 and B2 shocked regions. We report 
non-detections of NH 2 OH in both sources. We a perform non-LTE analysis of CH 3 OH 
observed in our CSO spectra to derive kinetic temperatures and densities in the shocked 
regions. Using these parameters, we derive upper limit column densities of NH 2 OH of 
< 1.4 X 10^^ cm“^ and < 1.5 x 10^^ cm“^ toward the B1 and B2 shocks, respectively, 
and upper limit relative abundances of Nj^ih^qh/Nh 2 < 1-4 x 10“® and < 1.5 x 10“®, 
respectively. 


1. Introduction 


Glycine (NH 2 CH 2 COOH), the simplest amino acid, has been identified in cometary (jElsila et al 


2009|) and meteoritic samples (jKvenvolden et al. Ill970l b however it is not yet known whether the 
species is formed in these solid bodies, i n the interstellar medium (ISM), or both. Indeed, a detection 


in th e ISM has so far proven ellusive (jSnvder et al. 1120051 : ICunningham et al. 1120071 : 1.1 ones et al 


20071 ). A lthough chemic al models predict the formation of glycine in the ISM, in low abundance 


(see, e.g. iGarrod (j2013l l and refs, therein), it has been difficult to observationally constrain them 
without a detection of the species. The relative likelihood of different possible formation routes 
can, however, be constrained by careful observation of the reactants used in the chemical models 
to form glycine. 

One pathway which has garnered significant interest is the formation of glycine through the 
reaction of hydroxylamine (NH 2 OH), or its protonated and ionized derivatives, with acetic acid 
(CH 3 COOH), a known interstellar molecule (jMehringer et al. 11199711. Ionization or protonation o f 


NH 2 OH under interstellar conditions should be efficient (|Angelelli et al. lll995l : lBoulet et al. Ill999l i 


and subsequent laboratory wo rk has demonstrated the formation of glyc ine from NH 2 OH and its 


ionized and protonated forms (|Blagoievic et al. 120031 : ISnow et al. 120071 1 . Although recent theo¬ 


retical work has suggested that these g as-phase routes forming glycine through these reactions are 


inefficient under interstellar conditions (|Barrientos et al. II2012I ) , condensed-phase surface-mediated 
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reactions remain a possibility, and observational constraints of the precursors are still desirable for 
the refinement of models. While acetic acid is readily-constrained, NH 2 OH has yet to be observed 
in the ISM. 

In recent years, laboratory work has shown that a number of pa t hways exist which result 


in efficient formation of NH 2 OH on grain surfaces. IZheng fc Kaiser (120101 ) demonstrated the 


formation of NH 2 OH in H 2 O-NH 3 ices after UV-irradiation. Formation by s uccessive hydrogenati on 


of NO has been shown to proceed efficiently, barrierlessly, and in h igh yield (ICongiu et al. 


Fedoseev et al. II2OI2I : IIoddoIo et al. II2014I : iMinissale et al. Il2014l i . Most recently, iHe et al. ( 201a) 


2012 a 




demonstrated the efficient formation of NH 2 OH via oxidation of NH 3 in ices with a low barrier to 
activation (of order 1000 K). 


Earlier modeling work bv iGarrod et al. ([200a) also suggested that NH 2 OH is formed in high 


abundance in interstellar ices, and is subsequently liberated into the gas-phase during the warm-up 
pe riod of emerging hot c ores and hot corinos in readily-detectable quantities. Observational efforts 
bv iPulliam et al. (j2012l ). however, failed to detect NH 2 OH toward a selection of such sources known 
to be rich in complex organic material (e.g. Sgr B2( N) and Orion-KL) est ablishing upper limits 


of NN HoOH/N ho < 10 ® — 10 Refinement of the iGarrod et al. (|2008l i model in subsequent 


work (Garrod 


20131 1. as well as limiting modeling presented in the aforementioned laboratory 


studies, resulted in p redicted gas-phase abundances in line with the upper limits established by 


Pulliam et al. ([2013). Yet, with the incl usion of the H + HNO —>■ HNOH pathway into the model. 


(| 2012 al ). condensed phase abundances of 


as described by the laboratory work of iGongiu et al. 

NH 2 OH again approach N]\fH 20 H/NH 2 ~ 10 “®, with gas-phase abundances of ^nh 20 h/^H 2 
10 -^. 


The iGarrod (|2013l l model largely assumes that the release of NH 2 OH into the gas-phase is 
a gradual process dominated by the warm-up of the hot core. Both the laboratory work and the 
models, however, predict NH 2 OH is initially formed in large abundance at very cold (< 20 K) 
temperatures and early in the evolution of these sources. Thus, the most optimistic source for a 
detection of gas-phase NH 2 OH is one where the reservoir of condensed-phase NH 2 OH formed at 
low temperatures is liberated en masse into the gas-phase prior to release by thermal mechanisms. 

Shocked regions displaying high degrees of molecular complexity likely represent this best-case 
scenario. In these regions, complex molecules are formed efficiently in ices at low temperatures, 
but are not otherwise liberated into the gas-phase except by thermal desorption at much greater 
temperatures. When these ices are subjected to shocks, however, the mantle is non-thermally 
ejected into the gas-phase, resulting in large abund a nces o f relatively cool {Trot < 100 K), complex 


organic material (see, e.g., iRequena- Torres et al. ([20061 )). One of the most prominent of these 


regions is the young protostellar outflow LI 157. Numerous recent studies report high degrees of 
molecular complexity arising from shocked regions within the outfl ow, which originates in cold, 
quiescent gas around the protostar (see, e.g.. lArce et al. ([2003) and ICodella et al. (l2015lB . 


Here, we present deep searches for NH 2 OH using the Caltech Submillimeter Observatory (CSO) 
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at A = 1.3 mm and the Combined Array for Research in Millimeter-wave Astronomy (CARMA) 
at A = 3 mm toward the L1157 outflow. We report non-detections of NH 2 OH in both searches. In 
order to derive upper limits to the column density of NH 2 OH, we use transitions of CH 3 OH observed 
with the CSO to constrain the kinetic temperature and density in the shocked gas using a radiative- 
transfer approach, and CARMA images of CH 3 OH to determine the size of the shocked gas. Finally, 
we estimate upper limits to the abundance of NH 2 OH and discuss possible implications. 


L1157 


L1157 is a dark cloud in Cepheus located a distance of ~250 pc (jLoonev et al. 1120071 1. and 
contains a prototypical shocked bipolar outflow from a Class 0 protostar. R has been the subject of 
great interest in the last twenty years, with numerous studies investigating the physical conditions 
within the source. Due to the variety of methods used for these studies, direct comparisons be¬ 
tween results is challenging. Nevertheless, an overall picture does come into view, and this general 
description will be sufficient for the discussion presented here. 


Originating in the cold, quiescent gas (T ~ 13 K, iBachiller et al. 


(j 19931 )1 surrounding the 

Class 0 protost a r L11 57-mm, the southern lobe of the accelerated outflow (T ~ 50 — 100 K, 


Bachiller et al. 


(jl993l ll has undergone two major shocks, referred to as L1157-B1 and L1157-B2. 
The B1 shock is younger and warmer (~2000 yr, T^in ~ 80 — 100 K) than the B2 shock (~4000 yr. 
Thin ~ 20 — 60 K), and many complex chemical species are observed in enhanced abundance towar d 
both shocks due to non-thermal desorption from grains ( Mendoza et al. l2014l : ICodella et ani201,5ll . 
While the absolute values of these physical parameters vary somewhat within the literature, this 
qualitative picture and enhancement in chemical abundance is consistently reported. 


3. CSO Observations 

The spectrum toward L1157 obtained with the CSO was collected over 8 nights in 2012 July, 
August, and September, and 7 nights in 2014 September and December. The telescope was pointed 
at the B1 shocked region at a(J2000) = 20^39™07®.7, (5(J2000) = 68°01'15".5 and the B2 shocked 
region at a(J2000) = 20^39™13®, (5(J2000) = 68°00'37" (see Figured]) and spectra were adjusted 
to a vlsr = 1-75 km s“^. A small subset of scans toward B1 were obtained at a second position 
offset by -4" in declination. The CSO 230/460 GHz sidecab double-sideband (DSB) heterodyne 
receiver was used in fair weather (r 22 o = 0.1 - 0.25) which resulted in typical system temperatures 
of 300 - 500 K. The backend was a fast Fourier-transform spectrometer which provided 1 GHz of 
DSB data at 122 kHz resolution (~0.2 km s“^ at 230 GHz). The total frequency coverage was 188 
- 193 GHz, 200 ~ 205 GHz, 237 - 243 GHz, and 249 - 255 GHz toward Bl, and 237 - 243 GHz 
and 249 ~ 255 GHz toward B 2 . The targeted transitions of NH 2 OH are given in Tabled! 

A chopping secondary mirror with a throw of 4' was used for ON-OFF calibration. Point- 
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ing was performed every ~2 hours and typically converged to within 1 ~ 2'h The raw data were 
calibrated using the standard chopper wheel calibration method resulting in intensities on the 
atmosphere-corrected temperature scale. These were then corrected to the main beam tem¬ 
perature scale, Tmb, where T^b = T^/rjmb- For these observations, the main beam efficiency was 
r]b = 0.70. The spectra were collected in DSB mode at a variety of IF settings to allow for a robust 
deconvolution. The CLASS package from the GILDAS suite of program^ was used for the data 
reduction. Spurious signals were removed from the spectra, which were then baseline subtracted 
using a polynomial fit. The standard CLASS deconvolution routine was used to generate single¬ 
sideband data. The spectra were then Hanning smoothed to a resolution of ~1.4 km s“^. The 
average full-width half-maximum linewidth in B1 was ~7.9 km s“^, and in B2 was ~5.2 km s“^. 


4. CARMA Observations 

A total of 89.3 hours of observations were conducted with the CARMA 15-element array in C- 
configuration (2013 May), D-configuration (2012 October, November), and E-configuration (2012 
August) at A = 3 mm. The phase center for these observations was a(J2000) = 20^39™07®.7, 
5(32000) = 68°01'll".5. The CARMA correlator was used in its 62 MHz bandwidth, 3-bit mode 
providing 255 channels across the band for a resolution of 243 kHz or 0.7 km s“^. The three 
targeted transitions of NH 2 OH are given in Table [H 

MWC349 and Neptune were used as primary flux calibrators; the passband calibrators were 
1635-^381, 2232-M17, 0102-^584, 1743-038, 2015-^372, and 3C84. The gain calibrator was 1927-^739. 
Data reduction was completed using the Miriad package. Data were flagged for antennas that were 
offline or malfunctioning during observation or when the phases in the calibrators had deviations 
greater than thirty degrees from a smooth trend line. Image cleaning was also performed using 
Miriad. The robust factor tended towards natural weighting with cell size set to 0.4" and image size 
of 1024 pixels. No channel averaging was applied for the maps shown; however, channel averaging 
to 1.4 km s“^ was used on windows targeting hydroxylamine transitions. Clean regions were drawn 
with the polygon tool around clear emission and cleaned to an average noise of 0.7 mJy beam“^. 
The restoring beam was typically ~3".4 x 3".2. 

The full set of observations toward L1157 CARMA included a host of other molecular species, 
and will be presented in a follow-up paper by Dollhopf et al. (2015). 


^Institut de Radioastronomie Millimetrique, Grenoble, France - http://www.iram.fr/IRAMFR/GILDAS 
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Table 1. 


Targeted NH 2 OH transitions, spectroscopic parameters, assumed collisional coefficients, 

and critical densities. 


Transition 

Frequency* 

(MHz) 

Eu 

(K) 

Sij fJ. 

(Debye^) 

(cm® 

7+ 

s-i) 

log(A) 

(log(s-i)) 

n, 

(cm 

t 

2i,2 

- ll,l 

100683.58(20) 

15.204 

0.520 

4.1 

X 

10 -“ 

-5.9079 

3 

X 

10"* 

2o,2 

- lo,l 

100748.23(20) 

7.2527 

0.694 

5.5 

X 

10-“ 

-5.7821 

3 

X 

lO"* 

2 i,i 

- ll,0 

100807.62(20) 

15.213 

0.520 

4.1 

X 

10-“ 

-5.9062 

3 

X 

10"* 

5i,5 

- 4i,4 

251677.3666(78) 

44.194 

1.665 

4.0 

X 

10 -“ 

-4.5515 

7 

X 

10® 

54,1 

“ 44,0 

251734.8061(82) 

163.51 

0.625 

2.4 

X 

10 -“ 

-4.9771 

4 

X 

10® 

54,2 

- 44,1 

251734.8061(82) 

163.51 

0.625 

2.4 

X 

10 -“ 

-4.9771 

4 

X 

10® 

53.3 

“ 43,2 

251780.3167(70) 

107.85 

1.110 

1.6 

X 

10 -“ 

-4.7270 

1 

X 

10® 

53.2 

- 43,1 

251780.3178(70) 

107.85 

1.110 

1.6 

X 

10 -“ 

-4.7270 

1 

X 

10® 

52,4 

“ 42,3 

251811.9913(70) 

68.081 

1.457 

2.0 

X 

10 -“ 

-4.6088 

1 

X 

10 ® 

52,3 

- 42,2 

251813.8609(70) 

68.081 

1.457 

2.0 

X 

10 -“ 

-4.6087 

1 

X 

10 ® 

5o,5 

“ 4o,4 

251838.4937(79) 

36.261 

1.735 

1.2 

X 

10-10 

-4.5329 

1 

X 

10 ® 

5i,4 

- 4i,3 

251987.1350(78) 

44.239 

1.665 

4.0 

X 

10 -“ 

-4.5498 

1 

X 

10 ® 


Transitions an d parameters acces sible at www.splata l oKue.n et. Original laboratory 


work reported by Tsunekawa ( 1972h and Mornio et al. 1 200d l . Catalogued at CDMS 


(Muller et al. 200 




^Values were taken from corresp onding transitions of ^-CHsOH obtained from the Leiden 


Atomic and Molecular Database iSchdier et al. 


tion. Original data reported bv iRabli fc Flower 
^ Given as n^r = Aj^ 


200^ and used without further modifica- 


( 20inli . 
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5. Non-LTE Modeling of CH 3 OH 


An accurate determination of NH 2 OH column density upper limits in the B1 and B2 shocks 
requires the kinetic temperature and density in these regions. To obtain these, we fit the observed 
CH 3 OH emission in the same spectral window (Table [2|) to a non-LTE model of the outflow and 
shock , based on a physical mo del derived from our our CARMA observations, using the RADEX 


code ( van der Tak et ah 1120071 ) 





Table 2. Observed CH 3 OH transitions, spectroscopic parameters, collisional coefficients, and 

critical densities. 


Transition^ 

Freqnency* 

(MHz) 

Eu 

(K) 

(Debye^) 

(cm^ 

7^ 

‘s-i) 

log(A) 

(log(s-i)) 

n * 

/ i-cr 

(cni“3) 




T-CH 3 OH 





2o,2 “ 

lo,l + + 

96741.375(5) 

6.9650 

1.617 

5.5 X 

10 -“ 

-5.4675 

6 X lO"! 

5i.5 - 

4i,4 + + 

239746.253 

49.059 

3.885 

4.3 X 

10 "“ 

-4.2468 

1 X lO’' 

5o,5 - 

4o,4 H" + 

241791.431 

34.817 

4.043 

1.2 X 

10-10 

-4.2184 

5 X 10® 

53,3 - 

43,2 H— h 

241832.91(20) 

84.618 

2.578 

1.6 X 

10-“ 

-4.4137 

2 X 10® 

53,2 “ 

43,1 - - 

241832.91(20) 

84.618 

2.578 

4.6 X 

10-“ 

-4.4137 

8 X 10® 

52,4 - 

42,3 - - 

241842.324 

72.530 

3.415 

1.3 X 

10-13 

-4.2915 

4 X 10® 

52,3 “ 

42,2 + + 

241887.704 

72.533 

3.415 

1.8 X 

10-13 

-4.2912 

3 X 10® 

5i,4 - 

4i,3 - - 

243915.826 

49.661 

3.885 

3.8 X 

10-13 

-4.2243 

2 X 10® 

llo,ll 

- 10i,io + + 

250506.980 

153.10 

10.63 

9.0 X 

10-13 

-4.0728 

9 X 10 '^ 

00 

w 

1 

82,6 - + 

251517.262 

133.36 

7.308 

8.1 X 

10-13 

-4.0990 

1 X 10® 

73,4 - 

72,5 - + 

251641.667 

114.79 

6.279 

1.4 X 

10-“ 

-4.1099 

6 X lO’' 

63,3 - 

62,4 —h 

251738.520 

98.546 

5.224 

6.7 X 

10-“ 

-4.1272 

1 X lO’' 

53,2 - 

52,3 - + 

251811.882 

84.618 

4.126 

1.8 X 

10-11 

-4.1567 

4 X 10® 

43,1 “ 

42,2 - + 

251866.579 

73.012 

2.954 

2.2 X 

10-11 

-4.2144 

3 X 10® 

53,3 - 

52,4 H— 

251890.901 

84.619 

4.125 

8.1 X 

10-“ 

-4.1564 

9 X 10® 

63,4 - 

62,5 + - 

251895.728 

98.545 

5.220 

2.6 X 

10-“ 

-4.1266 

3 X 10 '^ 

43,2 - 

42,3 + - 

251900.495 

73.012 

2.953 

4.5 X 

10-11 

-4.2143 

1 X 10® 

33,0 - 

82,1 - + 

251905.812 

63.727 

1.641 

3.1 X 

10-11 

-4.3603 

1 X 10® 

33,1 - 

82,2 + - 

251917.042 

63.727 

1.641 

1.3 X 

10-“ 

-4.3603 

3 X 10 '^ 

^3,5 “ 

72,6 + - 

251923.631 

114.79 

6.272 

8.5 X 

10-11 

-4.1089 

9 X 10® 

83,6 “ 

82,7 + - 

251984.702 

133.36 

7.293 

2.1 X 

10-11 

-4.0975 

4 X 10® 




T^-CHaOH 





2-1,2 

- 1-1,2 

96739.362(5) 

12.541 

1.213 

1.4 X 

10-13 

-5.5922 

2 X lO’' 

2o,2 - 

to ,2 

96744.550(5) 

20.090 

1.617 

9.5 X 

10-“ 

-5.4675 

4 X 10® 

2i,i “ 

li,o 

96755.511(5) 

28.011 

1.244 

3.7 X 

10-11 

-5.5811 

7 X 10"^ 

5o,5 - 

4o,4 

241700.219 

47.934 

4.040 

3.0 X 

10-11 

-4.2192 

2 X 10® 

5-1,5 

- 4-1,4 

241767.224 

40.391 

3.882 

7.9 X 

10-13 

-4.2362 

7 X lO’' 

53,2 - 

43,1 

241843.646 

82.531 

2.587 

1.5 X 

10-“ 

-4.4120 

3 X lO’' 

5i,4 - 

4i,3 

241879.073 

55.871 

3.980 

3.1 X 

10-“ 

-4.2248 

2 X lO’' 

5-2,4 

“ 4-2,3 

241904.152 

60.725 

3.399 

1.1 X 

10-13 

-4.2932 

5 X 10® 

52,3 “ 

42,2 

241904.645 

57.069 

3.356 

2.5 X 

10-“ 

-4.2987 

2 X 10 '^ 

2o,2 - 

1-1,1 

254015.340 

20.090 

0.499 

4.7 X 

10-13 

-4.7208 

4 X 10 '^ 
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Note. — Only transitions with modeled peak intensity above 5 mK are listed here. A complete 
modeled spectrum is available as Supplementary Information. 

“+” and refer to A+ and A— parity states, respectively. For iS-CHsOH, a negative value of 
Ka is used to differentiate between (+)i7i and (-)i ?2 states, which belong to the same E symmetry 
species. 


* Transi tions and parameter s accessible at www.splatalog ue.net. Original labo ratory work re¬ 
ported by Xu fc Lovas (1 1997 1 and references therein and by Muller et al. 1 2004 1. Catalogued at 


CDMS ( Muller et al. 2005h . Except where noted, uncertainties are 50 kHz. 


^Values obtained fr om the Leiden Atomic and Molecular Database i Schoier et al. 2005l l. Orig- 

(l2O10ll . 


inal data reported bv iRabli fc Flower 


^ Given as ricr = A/y 
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First, we determined the approximate spatial extent and location of the shocked gas from the 
CH 3 OH emission observed with CARMA. Figured] shows the CSO beam at the targeted positions 
overlaid on our CARMA observations of CH 3 OH and HNCC@ at A = 3 mm. The pointing position 
for B1 was chosen from the literature prior to our CARMA observations, and is roughly centered 
on the shock front, rather than the peak of the shock emission. The pointing position for B2 was 
chosen after we had acquired the CARMA observations, and is centered on the peak of the shocked 
HNCO emission. In both cases, the shock is nearly completely contained within the beam, and we 
estimate a reasonable equivalent size of ~ 12 " for the shock, and ~ 20 " for the second component, 
within our ~30" CSO beam. In Bl, the shock is not centered in the beam, and thus any emission 
will suffer more significant fall-off effects than in B2. 


We made the initial assumption that the two components in our fit represented the shocked gas 
and the outflow. Under these assumptions, and based on our CARMA observation as well as prior 
studies (|Bachiller et al. Ill995l b we constrained the fit of the CH 3 OH emission with the following 
assumptions: 


1. The kinetic temperature of the shock was no less than that of the outflow. 

2. The density of the shock was no less than that of the outflow. 

3. The column density of CH 3 OH was no less in the shock than in the outflow. 


Assumptions (1) and (2) are well-substantia ted in the literature, both tow ard L1157 ([Bachiller et al 


19951 1 and in shocked environments in general (Ivan Disheock &: Blake Ill998l l. and assumption (3) 
agrees with the previous findings of iBachiller et al. (Il995l ). Using these constraints, and the two- 
source component model described earlier, we performed a reduced-y^ analysis of a grid of models 
fo r the CSO observations , simultaneously fitting A- and U-CH 3 OH following the methods described 
Crockett et al. ( 20141 ). A thorough discussion of the uncertainties in the fitting method is pro- 


m 


vided in Appendix jA] 

Based on these uncertainties, and the assumptions described above, we find the parameters in 
Table [3] represent a best-fit to the data, with reduced-y^ values of 2.50 and 1.77 for Bl and B2, 
respectively. These values are likely due t o the simplicity of ou r model relative to the complexity 
of the source. Literature sources (see, e.g.. iLefloch et al. (j2012l )) suggest a structure with three or 
more components is more realistic, but such additional components would not be well-constrained 
in our case due to the limited number of observed CH 3 OH transitions. Additionally, there are 
likely quite large temperature and density gradients, especially in the newer, more compact Bl 
shock, which will contribute additional error to the fit which is not accounted for in our analysis. 
Finally, CH 3 OH excitation is sensitive to the far-IR radiation field present, which we have assumed 


^Details and analysis of the HNCO observations will be presented in Paper II by Dollhopf et al. (2015) 
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Fig. 1.— CARMA observations of CH 3 OH and HNCO toward L1157. Contours are an integrated 
Moment 0 map of the 2 _i ^2 — 2 o ,2 — lo,i++) and 2 o ,2 — lo,i transitions {E^ = 12.5 K, 7 K, 

and 20.1 K, respectively) of CH 3 OH at contour levels of 0.16, 0.32, 0.49, 0.66, and 0.82 Jy/beam. 
These are overlaid on an integrated Moment 0 map of the 4 o ,4 — 80,3 transition of HNCO spanning 
0 — 7 > 0.025 Jy/beam in a single color cycle. The HNCO transitions trace the warmer, compact 
shocks, while the CH 3 OH transitions reveal the overall structure of the colder, diffuse outflow. The 
location of the two targeted observations of B1 and B2 with the CSO are shown as dashed white 
lines, approximately equal to the CSO beam size (~30") at the observed frequencies. 
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as standard galactic background, any deviation from which would further impact the accuracy of 
the model and increased the ultimate 


Simulated non-LTE spectra of CH 3 OH from these results toward B1 and B2 are provided in Fig- 
ures[2]and[3l respectively, overlaid on observations and corrected for beam efficiency and a 12 " source 
size. The CH 3 OH column densities and H^;; densities found here are in relatively good agreement 
with previous observat ions of the source (|Bachiller et ah Ill995l : iBachiller &: Perez Gutierrez 1119971 : 


Sugimura et al. Il2010l b The derived values of Tkin for the w armer component, w hich we ascribe to 


the shocks, agrees well with previous measurements in B2 (|Lefioch et al. Il2012l b while the upper 
range of our derived value in B1 falls at the lower edge of previous measurements. 

The kinetic temperature of the second component {Tkin = 10 K) is significantly lower than 
previous measurements of both the shock and the outflow {Tkin ~ 60 — 120 K), indicating that it is 
likely not probing the outflow as we had originally presu med. Instead, toward Bl, this component 


is more likely related to the 53 component described by iLeffoch et al. (|2012l l as remnants from 


the gas in which the earlier B2 shock was formed, also offering a possible explanation for the 
component’s pres ence toward B2. The 93 component has been r eported to be both cold {Tkin = 23 


K) and extended iMendoza et al. 
with our fit. 


(|2014l i: iLefloch et al. (j2012l i. in qualitatively good agreement 


6 . Results 

Despite deep searches, we find no conclusive evidence for NH 2 OH emission at A = 3 mm 
or A = 1.3 mm in either the CARMA or the CSO data, respectively. To calculate appropriate 
upper limits for NH 2 OH in the both datasets, we assume that NH 2 OH, if present in the gas phase, 
originates in the warm, shocked gas traced by the warm, dense CH 3 OH component. The rationale 
for this assumption is discussed in ^ 


Table 3. Best-fit parameters for two components of CH 3 OH toward L1157-B1 and B2 
determined from non-LTE RADEX calculations. 


Parameter 

L1157-B1 

L1157-B2 

Component 1 

Component 2 

Component 1 

Component 2 

nH 2 (cm“®) 

3(1) X 10® 

3(1) X 10® 

6(2) X 10® 

6(2) X 10® 

Tkin (K) 

60(20) 

10(3) 

50(15) 

10(3) 

Nch^oh (cm“^) 

3(1) X 10^® 

5(2) X 10i'‘ 

3(1) X 10^® 

6(2) X 10^^ 


Note. — Errors given are from Equation I All and calculated to be ^32%. 
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Fig. 2.— Non-LTE simulation of CH 3 OH spectrum from RADEX fit in red, overlaid on CSO 
observations of the B1 shock in black. The 4 i ^4 — transition of A-CH 3 OH at 191810.5 MHz 

was not included as part of the RADEX fit. It is, however, reasonably well-reproduced by the fit 
of the higher-frequency lines, adding additional evidence for the robustness of the fit. 
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Fig. 3.— Non-LTE simulation of CH 3 OH spectrum from RADEX fit in red, overlaid on CSO 
observations of the B2 shock in black. 
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6.1. NHaOH in CSO Data 


Near the temperatures derived from our RADEX fit in the warm component toward B1 and 
B2, the strongest NH 2 OH transition is the 60,5 - 4 o ,4 at 251838 MHz. As noted in Tabled] the 
critical density {ricr) for this transition is 10® cm“®. This value is based on the assumption that 
the collisional coefficients for NH 2 OH are similar to those of A-CH 3 OH. While not exact, the 
similar mass, dipole moments, molecular size, and energy level structures make this a reasonable 
approximation within the context of the following discussion. 


Given the densities derived from our RADEX ht (3 - 6 xlO® cm“^), we therefore assume that, 
at least for this transition, LTE is a reasonable approximation for determining NH 2 OH upper limits 
and thus T^x = Tkin- Under these conditions, we derive Icr upper limit column densities for NH 2 OH 
of < 1.4 X 10^“^ cm“^ and < 1.0 x 10^^ cm“^ in B1 and B2, respectively. 


Lefloch et al. ([2013) derive a CO column density of 0.9 x 10 


17 2 


in the §2 compo¬ 
nent of the B1 shock, arising from the shocked gas and covering the entire region, assuming a 
20'' source size. Assuming the CO is homogeneously distributed over the region, this gives a 
source-averaged C O column of 3.2 x 10^ ® cm ~^ for the 12" s ource size used here. Taking Nh 2 = 
3 X 10^(Nco) (c-f. Bolatto et al. (2013) and Lefloch et al. ( 2012 ')!. this results in an H 2 column 
density o f ^ 10 ^^ cm~^. A similar argument u sing the CO column density of 5.2 x 10^® cm~^ for B 2 
found bv iBachiller fc Perez Gutierrez ( 19971 ) also results in an H 2 column density of 10 ^^ cm 
This gives upper limits to the fractional abundances of NH 2 OH of N]\[H 20 h/]^H 2 1-4 x 10“^ 

and < 1.0 X 10“^ for B1 and B2, respectively. 


6.2. NH 2 OH in CARMA Data 

To calculate appropriate upper limits for NH 2 OH in the CARMA data, we again assume that 
NH 2 OH, if present in the gas phase, originates in the warm, shocked gas traced by the warm, 
dense CH 3 OH component. The spectra extracted from a 12" beam centered on these locations are 
shown in Figured] While there is an unidentified feature at the frequency of the 2 i ^2 ~ li,i NH 2 OH 
transition, if this were truly NH 2 OH emission, the other transitions in the observed window would 
be of equal or greater intensity. 

We measure the RMS to be 1.6 mK and 2.9 mK in the B1 and B2 spectra. At the derived 
temperatures, the 2 o ,2 — lo,i transitions at 100748 MHz is predicted to be strongest. As both B1 
and B2 have densities higher than Ucr for this transition (see Table [T]), we again derive an upper 
limit to the column density assuming LTE conditions. 

Following the same procedure as for the CSO data, we find a Icr upper limit of NH 2 OH of 
< 1.4 X 10^® cm“^ for B1 and < 1.5 x 10^® cm“^ for B2. This gives upper limits to the fractional 
abundances of NH 2 OH of Nj^h^qh/Nh 2 =< l-4x 10“® and < 1.5x10“® for B1 and B2, respectively. 
All results are summarized in Table d] 














- 16 - 



Fig. 4.— Spectra toward L1157-B1 (bottom) and B2 (top) extracted from 12" beams centered 
on the peaks of the shocked CH 3 OH emission and smoothed to a velocity resolution of ~1.4 km 
s“^. Red lines indicate the frequencies of NH 2 OH transitions in this frequency region. A U-line is 
present toward the B1 shock at the frequency of the 2 i ^2 ~ li,i NH 2 OH transition. 


Table 4. Derived upper limits to NH 2 OH column density and relative abundance in L1157-B1 

and L1157-B2 from the CSO and CARMA data. 



CSO A 

= 1 mm 

CARMA A 

= 3 mm 


B1 

B2 

B1 

B2 

NNH 2 OH (cm 
Nnh^oh /Nh2 

< 1.4 X 10^'* 

< 1.4 X 10"^ 

< 1.0 X 10^'‘ 

< 1.0 X 10"^ 

< 1.4 X 10^^ 

< 1.4 X 10"® 

< 1.5 X 10“ 

< 1.5 X 10"® 
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Discussion and Conclusions 


In ^ we assumed that NH 2 OH, if present, w ould arise predominant ly the warm, shocked 
regions B1 or B2, rather than the molecular outflow. IZheng fc Kaiser ([20101) show that NH 2 OH is 
thermally liberated from their laboratory samples between 160 - 180 K, whe reas for realistic inter - 
stellar ices and conditions, temperatures above ~110 K are likely sufficient (ICollings et al. II2004I ). 


Both our non-LTE fits, and the literature values discussed earlier have shown that temperatures in 
the targeted regions are below these thresholds, and thus a thermal mechanism for the desorption 
of NH 2 OH from grain surfaces in L1157 is unlikely. Thus, non-thermal desorption in the shocks 
should be the dominant mechanism for NH 2 OH liberation, and any NH 2 OH emission should trace 
these shocked regions. 

The upper limit column d e nsitie s established with the CARMA data are equivalent to those 
established by Pulliam et al. ( 2012 1 of < (0.9 — 8) x 10^^ cm~^. The limits established by 


Pulliam et al. (|2012l l. however, are beam-averaged column densities, whereas the CARMA mea¬ 


surements presented here provide a more well-constrained source size, and thus a more robust upper 
limit. Nevertheless, t he upper limits t o rela tive abundance found here are significantly higher than 


those determined by iPulliam et al. (120121 1 of ~10 due to the low H 2 column in the region. 


These low columns are common in ot her chemically-complex shock ed regions as well, such as a host 


of Galactic Center clouds studied bv iReouena- Torres et al. (j2006l l (N wr, = (2 — 68l x 10^^ cm 


and other shocked outflows like BHR 71 ~ H) ^ 10^^ cm ^: lGarav et al. (1998])) 


Yet, these regions likely represent the best-case scenario for a gas-phase detection of NH 2 OH. 
Such a detection is critical to accurately constrain the application of laboratory results, which show 
NH 2 OH is a significant player in grain-surface nitrogen chemistry, to chemical models. The deep 
search presented here stretched the capabilities of the CSO and CARMA to the limit; searches 
deep enough to be sensitive to the relative abundances predicte d by laboratory work and modeling 
are simply time-prohibitive on these types of instruments. The lGarrod (120131 1 model, without the 
addition of the H -I- HNO —> HNOH pathway, predicts peak condensed-phase NH 2 OH abundances of 
~(5 — 10) X 10“®, which, if liberated en masse in a shock, is a detectable population, but likely only 
with the sensitivity and spatial resolution of ALMA. The upper limits established here, however, 
are several orders of magnitude lower than both the condensed-phase and gas-phase populations 
predicted by the augmented model which includes this hydrogenation pathway. Thus, our results 
demonstrate that further efforts are needed both in modeling and in the laboratory, to identify and 
fully-constrain both the formation and the destruction pathways for NH 2 OH. 


BAM gratefully acknowledges funding by an NSF Graduate Research Fellowship during initial 
portions of this work. The authors thank the anonymous referee for comments which increased 
the quality of this work. We thank S. Radford for support of the GSO observations. The National 
Radio Astronomy Observatory is a facility of the National Science Foundation operated under 
cooperative agreement by Associated Universities, Inc. Portions of this material is based upon 



































- 18 - 


work at the Caltech Submillimeter Observatory, which was operated by the California Institute 
of Technology under cooperative agreement with the National Science Foundation (AST-0838261). 
Support for CARMA construction was derived from the Gordon and Betty Moore Foundation, the 
Kenneth T. and Eileen L. Norris Foundation, the James S. McDonnell Foundation, the Associates 
of the California Institute of Technology, the University of Chicago, the states of California, Illinois, 
and Maryland, and the National Science Foundation. 


A. Error Analysis 


The uncert a inty i n the observations is given as Equation lAll modified from Equation A1 in 


Crockett et al. ( 2014l l. 


CTtot — \J Cr^ms + <^cal + + ^bf 

arms is the uncertainty due to the RMS noise level of the spectra. This value is consistently 
~0.0058 K across the CH 3 OH transitions in Bl, and ~0.0078 K in B2. 

a^ai is the uncertainty in the absolute flux calibration of the observations, which is ~30%. 

apt is the uncertainty in the pointing accuracy, which is taken to be 2 " based on the average 
convergence of pointing observations throughout the observing period. 

at,f is the uncertainty in the beam-hlling factor. The size of the telescope beam varies by 
~4" across the band, is well-described, and is calculated independently at each transition frequency, 
thus we assume no contribution to a^f from the telescope beam. However, the sizes of the two source 
components are not a varied parameter, and are only modestly well-described by our CARMA 
observations. Further, in the case of Bl, the shock is positioned off-center of the beam, making the 
emission more sensitive to fall-off effects that are not accounted for explicitly in the calculations. 
We therefore assume a factor of two uncertainty in source size to take both issues into consideration, 
which is reflected in a^f. 

These result in net uncertainties of ~32% in the absolute intensity of the observed transitions, 
with the largest contribution from the absolute flux calibration, values outside of the extremes 
represented by these uncertainties more than double those of the best-fit parameters, so we take 
~32% to be a reasonable estimate of the overall uncertainty in n// 2 ) NCH 3 OH, and T^in- 
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